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Abstract 

We study a system of N interacting particles on Z. The stochastic dynamics consists of 
two components: a free motion of each particle (independent random walks) and a pair-wise 
interaction between particles. The interaction belongs to the class of mean-field interactions 
and models a rollback synchronization in asynchronous networks of processors for a distributed 
simulation. First of all we study an empirical measure generated by the particle configuration 
on R. We prove that if space, time and a parameter of the interaction are appropriately scaled 
(hydrodynamical scale), then the empirical measure converges weakly to a deterministic limit 
as N goes to infinity. The limit process is defined as a weak solution of some partial differential 
equation. We also study the long time evolution of the particle system with fixed number of 
particles. The Markov chain formed by individual positions of the particles is not ergodic. 
Nevertheless it is possible to introduce relative coordinates and prove that the new Markov 
chain is ergodic while the system as a whole moves with an asymptotically constant mean 
speed which differs from the mean drift of the free particle motion. 

MSC 2000: 60K35, 60J27, 60F99. 

1 Introduction 

We study an interacting particle system which models a set of processors performing parallel sim- 
ulations. The system can be described as follows. Consider N > 2 particles moving in Z. Let 
Xi(t) be the position at time t of the i— th particle, 1 < i < N. Each particle has three clocks. 
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The first, the second and the third clock, attached to the i— th particle, ring at the moments of time 
given by a mutually independent Poisson processes ilj a , H^p and Ilj jAiJV with intensities a, (3 and 
/xjv correspondingly. These triples of Poisson processes for different indexes are also independent. 
Consider a particle with index i. If the first attached clock rings, then the particle jumps to the 
nearest right site: Xi — ► %i + 1, if the second attached clock rings, then the particle jumps to the 
nearest left site: Xi — > Xi — 1. At moments when the third attached clock rings a particle with index 
j is chosen with probability 1/N and if X{ > Xj, then the i— th particle is relocated: X{ — > Xj. It is 
supposed that all these changes occur immediately. 

The type of the interaction between the particles is motivated by studying of probabilistic mod- 
els in the theory of parallel simulations in computer science ( ifToll . Il2"2l l2"31l and O HB). The 
main peculiarity of the models is that a group of processors performing a large-scale simulation 
is considered and each processor does a specific part of the task. The processors share data do- 
ing simulations therefore their activity must be synchronized. In practice, this synchronization 
is achieved by applying a so-called rollback procedure which is based on a massive message ex- 
change between different processors (see [2, Sect. 1.4 and Ch. 8]). One says that Xi(t) is a local 
time of the i-th processor while t is a real (absolute) time. If we interpret the variable Xi(t) as 
an amount of job done by the processor i till the time moment t, then the interaction described 
above imitates this synchronization procedure. Note that from a point of view of general stochastic 
particle systems the interaction between the particles is essentially non-local. 

We are interested in the analysis of asymptotic behaviour of this particle system. First of all 
we consider the situation as the number of particles goes to infinity. For every finite A" and t we 
can define an empirical measure generated by the particle configuration. It is a point measure with 
atoms at integer points. An atom at a point k equals to a proportion of particles with coordinate 
k at time t. It is convenient in our case to consider an empirical tail function corresponding to 

1 N 

the measure. It means that we consider £,x,N{t) = -77 z2 1 (xi(t) > x) the proportion of particles 

having coordinates not less than x G R. The problem is to find an appropriate time scale t N 
and a sequence of interaction parameters fi^ to obtain a non-trivial limit dynamics of the process 
^N,[xN](tN) as A" — > 00. The cases a 7^ (3 and a = (3 require different scaling of time and the 
interaction constant ji^. We prove that there exist non-trivial limit deterministic processes in both 
cases as A" goes to infinity if we rescale time and the interaction constant as t N = tN, [i N = fj,/N 
in the first case and as t N = tN 2 , fi N = [i/N 2 in the second case respectively. The processes 
are defined as weak solutions of some partial differential equations (PDE). It should be noted that 
the PDE relating to the zero drift situation is a famous Kolmogorov-Petrowski-Piscounov-equation 
(KPP-equation, 0). This result was announced in ifP^ll . 

Another issue we address in the paper is studying of the long time evolution of the particle 
system with fixed number of particles. It is easy to see that the Markov chain x(t) = {xi(t), i = 
1, . . . , A^}, t > 0, is not ergodic. Nevertheless the particle system possesses some relative stabil- 
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ity. We introduce new coordinates yi(t) = Xi(t) — min, Xj(t), i = 1, . . . , N, and prove that the 
countable Markov chain y(t) = {yi{t), i = 1, . . . , iV}, t > 0, is ergodic and converges expo- 
nentially fast to its stationary distribution. Therefore the system of stochastic interacting particles 
possesses some relative stability. We show also that the center of mass of the system moves with 
an asymptotically constant speed. It appears that due to the interaction between the particles this 
speed differs from the mean drift of the free particle motion. 

It should be noted that the choice of the interaction may vary depending on a situation. Various 
modifications of the model can be considered and similar results can be obtained using the same 
methods. We have chosen the described model just for the sake of concreteness. 

Probabilistic models for parallel computation considered before by other authors. The pa- 
per [ 16 1 deals with a model consisting of two interacting processors (N = 2). It contains a rigor- 
ous study of the long-time behavior of the system and formulae for some performance character- 
istics. Unfortunately, there are not too many mathematical results about multi-processor models 
( lfT2l IT3l fT4l l22l 123 1). Usually mathematical components of these papers have a form of prepara- 
tory considerations before some large numerical simulation. The paper J6) is of special interest 
because it rigorously studies a behavior of some model of parallel computation with N processor 
units in the limit N — > 00. A stochastic dynamics of |6| is different from the dynamics studied in 
the present paper and main results of [6| concern a so-called thermodynamical limit. The authors 
prove that in the limit the evolution of the system can be described by some integro-differential 
equation. In the present study we propose a model which dynamics is easy from the point of view 
of numerical simulations and, at the same time, provides us with a new probabilistic interpretation 
of some important PDEs including the classical KPP-equation. 

The paper is organised as follows. We formally define the particle system, introduce some 
notation and formulate the main results in Sectional Sections and HI contain the proofs of the 
main results. In Section|5]we discuss solutions of the limiting equations. 

Acknowledgments. We are thankful to Dr. T. Voznesenskaya (Faculty of Computational Math- 
ematics and Cybernetics, Moscow University) who first introduced us to stochastic algorithms for 
parallel computations. The authors would like to thank Prof. V. Bogachev (Faculty of Mechanics 
and Mathematics, Moscow University) for the helpful discussions on convergence of measures on 
general topological vector spaces and for the suggested references. We are also grateful to Prof. 
V. Malyshev for his warm encouragement and valuable comments on the present manuscript. 
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2 The model and main results 

Formally, the process x(t) = {xi(t),i = 1, . . . , TV}, describing positions of the particles, is a 
continuous time countable Markov chain taking values in Z N and having the following generator 

N 

G N g{x) = (« (g (x + e™) - g(xj) + (3 (g (x - e\ N) ) - g(x))) + 
i=i 

N 

i=l j^i 

where x = (xi, . . . ,x N ) G Z^, <? : Z w — > R is a bounded function, e\ N ^ is a iV-dimensional 
vector with all zero components except %— th which equals to 1, /{ Xi>Xj .} is an indicator of the set 

{Xi > Xj}. 

Define 

1 N 

&v,fc(*) = y^^(t)>fc}, /c G Z. (2) 
i=i 

The process £/vW = {6v,fc(^) ; A: G Z} is a Markov one with a state space ifjv the set of all non- 
negative and nonincreasing sequences z = {z k , k G Z} such that z fe G {//-AT, / = 0, 1, . . . , N} for 
any k G Z and 

lim Zfc = 1, lim ^ = 0. 

fc— +— oo fc-++oo 

The generator of the process £at(0 is given by the following formula 

Ljv/(z) = N^dfiz + e fc /iV) - f(z))a(z k ^ - z k ) + (/(z - e k /N) - f{z))(5{z k - z k+1 )) 
k 

(3) 

+ N hn - (e J+ i + . . . + e fc )/JV) - f(z))(z k - z k+1 )(zi - z l+1 ) , 

Kk 

where e i: i G Z is an infinite dimensional vector with all zero components except %— th which 
equals to 1, / : H N — > R is a bounded function. 

2.1 Hydrodynamical behavior of the particle system 

Denote C,N,x{t) = ^N,[Nx](t), x G R. The process Cjv(^) takes values in H = H(R) the set of all 
non-negative right continuous with left limits nonincreasing functions having the following limits 

lim ip(x) = 1, lim ip(x) = 0. 
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Denote by S(R) the Schwartz space of infinitely differentiable functions such that for all m,n E 
Z+ 

||/|| m , n = sup|x m /W(x)| <oo. 

Recall that S(R) equipped with a natural topology given by seminorms || ■ || m) „ is a Frechet 
space (CHI). 

Define for every h E H a functional 

(h,f) = J h(x)f(x)dx, /6 5(H), 

R 

on the Schwartz space 5(R). The following bound yields that for each h E H (h, ■) is a continuous 
linear functional on 5(R) 

|(/>,/)|< / |/(x)|i^^<7r(||/|U + ||x 2 /||oc)=7r(||/|| 0> o+||x 2 /|| 2i o), 
7 1 + ar 

R 

where || • is the supremum norm. Thus the set of functions H(R) is naturally embedded into 
the space of all continuous linear functionals on 5(R), namely into the space S'(R) of tempered 
distributions. We will interpret (jv(£) as a stochastic process taking its values in the space S'(R). 

There are two reasons for embedding H (R) into S"(R) and considering the S"(R) -valued pro- 
cesses. The first reason is that due to some nice topological properties of S"(R) we can use in 
Section E] many powerful results from the theory of weak convergence of probability distributions 
on topological vector fields. And, secondly, the choice of S"(R) as a state space is convenient from 
the point of view of possible future study of stochastic fluctuation fields around the deterministic 
limits obtained in our main theorem |2~T1 

In the sequel we mainly deal with the strong topology (s.t.) on S"(R) (see Section |A~TT) . 
From now we fix some T > and consider ( N as a random element in a Skorokhod space 
D([0, T], 5"(R)) of all mappings of [0, T] to (S"(R), s.t.) that are right continuous and have left- 
hand limits in the strong topology on £'(R). Note that (S'(K),s.t.) is not a metrisable topological 
space therefore it is not evident how to define the Skorokhod topology on the space D( [0, T] , S'(R) ) . 
To do this we follow [TT31 and refer to Section lATTl 

Now we are able to consider probability distributions of the processes ((N(tN a ),t E [0, T}), 
a = 1, 2, as probability measures on a measurable space (-D([0, T], S"(R)), Bd([o,t],s'(r))) where 
Bd([o,t],s>(r)) is a corresponding Borel cr-algebra. It was proved in |7 1 that <6 d([ o,t],s'(r)) = C d([0 ,t],s'(r)), 
where Cd([o,t],S'(r)) is a cr-algebra of cylindrical subsets. 

Consider two following Cauchy problems 

u t (t, x) = -Xu x (t, x) + ^{u 2 (t, x) - u(t, x)) , (4) 
u(0, x) = ip(x) 
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and 

ut(t, x) = ju xx (t, x) + fi(u 2 (t, x) - u(t, x)) , (5) 
u(Q, x) = ip(x) 

where u t , u x and u xx are first and second derivatives of u with respect to t and x. Notice that the 
equation © is a particular case of the famous Kolmogorov-Petrowski-Piscounov-equation (KPP- 
equation, |9|). We will deal with weak solutions of the equations © and © in the sense of 
Definition EUl 

Fix T > and denote by C™ T = C£°([0, T] x R) the space of infinitely differentiable functions 
with finite support and equal to zero for t = T. 

Definition 2.1 (i) The bounded measurable function u(t, x) is called a weak (or generalized) so- 
lution of the Cauchy problem (0) in the region [0, T] x R, if the following integral equation 
holds for any function f G Cq° t 

T 

u(t, x)(f t (t, x) + \f x (t, x)) + fxu(t, x)(l — u(t, x))f(t, x))dxdt 



+ / u(0,x)f(0,x)dx = 




R 



(ii) The bounded measurable function u(t,x) is called a weak (or generalized) solution of the 
Cauchy problem (0) in the region [0, T] x R, if the following integral equation holds for any 
function f G C^ T 

T 

u(t,x)(f t (t,x) +jfxx(t,x)) +fm(t,x)(l - u(t,x))f(t,x))dxdt 



+ / u(0,x)f{0,x)dx = 




R 



R 

In Subsection l3.5l we will show that the both of Cauchy problems @ and © have unique weak 
solutions in the sense of Definition ^. II Here we want just to mention that this problem is not trivial. 
Indeed, the equation © is an example of a quasilinear first order partial differential equation. It is 
known that in a general case such type of equations might have more than one weak solution and 
it is only possible to guarantee uniqueness of the solution which satisfies to the so-called entropy 
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condition. The most general form of this condition was introduced by Kruzhkov in ifTUl . where 
he also proved his famous uniqueness theorem. Fortunately, in our particular case of the equation 
© the situation is quite simple due to simplicity of characteristics, they are given by the straight 
lines x(t) = Xt + C, do not intersect with each other and do not produce the shocks. Detailed 
discussions of the problem of uniqueness for equations @ and © are presented in Sub section Q3] 
The first theorem we are formulating describes the evolution of the system at the hydrodynam- 
ical scale. 

Theorem 2.1 Assume that an initial particle configuration £tv(0) = {^at,a, (0), k G Z} is such that 
for any function f G S(R) 

i im TfJ2Uk(0)f(k/N) = [ i>(x)f(x)dx, (6) 

N^oo iV z — ' J 
k R 

where ip G H(R). 

(i) If a — (3 = X 7^ and /jn = ju/iV, then the sequence {Qn-\}n=2 of probability distributions of 

random processes {(iy(tN), t G [0, T]}^ =2 converges weakly as N — » oo to the probability 
measure on D([0, T], S'(R)) supported by a trajectory u(t, x), which is a unique weak 
solution of the equation @ with the initial condition u(0, x) = ip(x) and as a function of x 
u{t,-) G H(R)Joranyt> 0. 

(ii) If a = (3 = 7 > 0, /i^ = n/N 2 , then the sequence {Qn\}n=2 of probability distributions of 

random processes {Qn if N 2 ), t G [0,T]}^ =2 converges weakly as N — > oo to the probability 
measure on D([0, T], S'(R)) supported by a trajectory u(t, x), which is a unique weak 
solution of the equation © with the initial condition u(0, x) = ip(x) and as a function of x 
u(t,-) G H{R)Joranyt> 0. 

2.2 Long time behavior of the particle system with the fixed number of par- 
ticles 

The number of particles is fixed in this section. Consider the following stochastic process y[t) = 
(yi(t),...,y N (t)), where 

yi(t) = Xi(t) - mmXj(t). 

j 

Note that x^ — xi = yu~ Vi for any pair k, I. 

It is easy to see that y(t) is a continuous time Markov chain on the state space 

r = |Jr fe cz^ 

k 
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where F k := {(z x , . . . , z k -i, 0, z k+1 , ...,z N ): z 3 e Z+}. 

Theorem 2.2 The Markov chain (y(t),t > 0) is ergodic and converges exponentially fast to its 
stationary distribution 

£ \P(y(t) = y)- 7t(y)\ < C ie x P (-C 2 t) 

uniformly in initial distributions ofy(0). 

3 Proof of Theorem O 
3.1 Plan of the proof 

The proof of the convergence uses the next well-known general idea (see, for example, ll2Tl § 5]). 
Let {a n } be a sequence in some Hausdorff topological space and assume that {a n } satisfies to the 
following two properties: (a) for any subsequence of {a n } there is a converging subsequence (this 
property is called a sequential compactness); (b) {a n } contains at most one limit point. Then the 
sequence {a n } has a limit. 

In our situation the role of {a n } is played by the sequences {Q^\}n=2 an d {Qn\}n=2- O ur 
proof consists of the following steps. 

Step 1. We fix an arbitrary T > and prove that the sequences of probability measures 

{Q < n\}n=2 an d {Qn-\}n=2 are tight- We use the Mitoma theorem ( lfT51 ) and apply martingale 
techniques widely used in the theory of hydrodynamical limits of interacting particle systems 

(f&ED). 

It is important to note that if a topological space V is not metrisable then, generally speaking, 
the tightness of a family of distributions on V does not imply a sequential compactness (see, for 
example, [3, V. 2, § 8.6]). So, in general, the above property (a) does not follow directly from the 
step 1. But in our concrete case V = D([0, T], S"(R)) we can proceed as follows. It was shown 
in that any compact subset of D([0, T], S"(R)) is metrisable. Due to this property we can apply 
the theorem from lETl Th. 2, § 5] which states that (under assumption of metrisability of compact 
subsets) the tightness of a family of measures implies its sequential compactness. All this justifies 
the next step. 

Step 2. We show that a measure that is a limit of some subsequence of the sequence {Q^ \}n=2 

(or {Q ( n\}n=2^> i s supported by the weak solutions of the partial differential equation © (or, 
correspondingly, ©). Then we note that each of the equations © and © has a unique weak 
solution (Subsection l3.5l) . This gives the above property (b). 
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3.2 Technical lemmas 

We start with some bounds which will be used throughout the proof. Denote 

k 

for z G Fat and / G 5(R). 

Lemma 3.1 (i) Ifaj^/3 and hn = fJ>/N, then for any z G if at 



C 
iV' 



< (V) 



iV - 2R f {z)L N R f (z)) = O (J^j . (8) 

(ii) If a = j3 and fi^ = fi/N 2 , then for any z G H N 



\L N Rf(z)\ < ^, (9) 



and 



iV 2 (L*J2£(z) - 2i2 / (z)L Ar it! / (z)) = O Qj) . (10) 
Zn Z?of/z cases C = C(f, a, (3, //). 

Proof of Lemma EU We will prove the bounds Q and (HI), the other ones can be proved simi- 
larly. We start with the bound Q. Using the equations 

R f (z + e k /N)-R f (z)= /U ' /A) 



R f (z - e k /N) - R f (z) 



N 2 ' 



AT2 

we get that for every z G .f/jv 



L N R f (z) =1 £ **(/?/((* - l)/iV) - (a + 0)/(*/JV) + a/((* + 1)/JV) 



if 

AT2 



' ^2f{k/N)z k {l-z k ) 
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For any function f E S (R) consider its upper Darboux sum 

U+(f) = ^-V max f(y). (11) 



Since U^(f) — > J f(x) dx as iV — > oo, the sequence {^jv(/)}jv=i * s bounded in AT for any fixed 
/. We have uniformly in z E H N 

\L N R f (z)\ < l(j2Hf((k + l)/N)-f(k/N))+(3(f(k/N)-f((k-l)/N))\ 



k 



< ^(l«-/9W(l/r|)+/^it(l/D), 

where / x = df(x)/dx. So the bound is proved. 

Let us prove the bound ©. Note that Ln = L$ + L$, where 

I$f(z) = Nj2(f(z + e k /N) - f(z))a(z k ^ - z k ) + 

k 

N J2(f(z - e k /N) - f{z))(3{z k - z k+1 ) 

k 

and 

Using the equations 
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one can obtain that for any z e H N 

LfR%z) = aNj2( R % z + tk/N) - R%z)){z k ^ - z k ) 

k 

+ /3NJ2( R2 A Z - e k/N) - R){z)){z k - z k+1 ) 



k 

: 2 M*)% E f(k/N)(z k ^ - Zk ) + ^Y, / 2 ( - **) 

k k 

2R A z ^ E /(*/#)(** - + ^ E fwwzk - 
it k 



2R f (z)L i ° ) R f (z) + 0[^). (12) 



1 

Direct calculation gives that for any function g k ,j(z) = z k Zj, k < j, on Hn 

L SvViO) = -jj( z kZj(l - z k ) + z k Zj(l - Zj)) + j^Zj(l - z k ). (13) 



Using this formula we get that for any z G H 



N 



L^Rj(z) = -2R f (z)^ Yl - **) 

k 

+ % E /(w/oviv)^(i -^) + |jE / 2 (w^(! - 

k<j k 



Summing the formulas (fT2l and (fl4ll we obtain 

Ljv^) = 2Rf(z)L N R f (z) + O (J-^j . (15) 

Lemma l3~T1 is proved. 



3.3 Tightness 

We will make all considerations for the part (i) (a ^ (3 and [i = fi/N) of the theorem. All reason- 
ings and conclusions are valid for the part (ii) (a = (3 and [i = n/N 2 ) of the theorem with evident 



changes. We will denote by x the probability distribution on the path space D[[0,T\, H 



A' J 



corresponding to the process ^{t) and by \ the expectation with respect to this measure. 
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Theorem 4.1 in lfT5ll (see also Section lA~2l yields that tightness of the sequence of {Q% a}n=2 



will be proved if we prove the same for a sequence of distributions of one-dimensional projection 
{(( N (tN), f),te [0, T]}^ =2 for every / E S(R). So fix / E S(R) and consider the sequence of 
distributions of the processes (^(tN), f), t E [0, T\. Note that the probability distribution of a 
process (^(tN), f) is a probability measure on D([0, T], R) the Skorokhod space of real-valued 
functions. 

By definition of the process (n^N) we have that 

(k+l)/N 

(( N {tN), f) = J2Uk(tN) J f(x)dx. 

k k/N 

It is easy to see that 

(( N (tN)J)=R f (Z N (tN)) + Mt)i 

where the random process 0at (if:) is bounded 1 0Ar(t) | < C(f)/N for any t > and sufficiently large 
N. Therefore it suffices to prove tightness of the sequence of distributions of random processes 
{R f (£ N (tN)),tE [0,7TO =2 . 
Introduce two random processes 

t 

W f)N (t) = Rf(£ N (tN)) - fl/(6r(0)) -N J L N R f (£ N (sN))ds, (16) 



and 



where 



V f , N (t) = (W f , N (t)Y - J Z f , N (s)ds, 

o 



Z LN (s) = N (L N R}(£ N (sN)) - 2R f (£ N (sN))L N R f (£ N (sN))) . (17) 

It is well known (Theorem 2.6.3 in [4| or Lemma Al.5.1 in ||8j|) that the processes Wf^if) and 
Vf t N(t) are martingales. 
The bound Q yields that 



N [ L N R f {£ N {sN))ds 



< C9 (a.s. 



for any time moment r. Thus the sequence of probability distributions of random processes 
{N f Q * L N R f (£ N (sN))ds, t E [0, T]}^ =2 is tight by Theorems CO and from Appendix. 
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The bound © yields that 



E^l(W LN (r + 0)-W fiN (r)) 2 = E 




Z ftN (s)ds < 



ce 



(18) 



for any stopping time r > since Vf^(s) is martingale. Using this estimate and Chebyshev 
inequality we obtain that the sequence of probability distributions of martingales {Wf^(t), t £ 
[0,T]}^ =2 is also tight by Theorem IA.31 Thus the sequence of probability distributions of the 
processes {Rf(^ N (tN)), t e [0, T]}^ =2 is tight by the equation (fToT) and the assumption © and, 
hence, the sequence of probability measures {Q^\}n=2 * s tig nt by Theorem 4.1 in lTT3Tl . 

3.4 Characterization of a limit point 

We are going to show now that there is a unique limit point of the sequence {Q ( n\} < n=2 an d this 
limit point is supported by trajectories which are weak solutions of the partial differential equation 
© in the sense of Definition l2.11 

Let f(s,x) G Cqj, and denote 



R,(t,MtN)) 



N ^ 

k 



Uk(tN)f(t,k/N), 



Define as before two random processes 




o 



and 




o 



where 



Z' LN (s) = N (L N R 2 f (s,£ N (sN)) - 2R f (s,Z N (sN))L N R f (s,Z N (sN))) . 
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By Lemma A 1.5.1 in [8| the processes Wj N (t) and Vj N (t) are martingales. It is easy to see that 



W' m {t) = (C(tN), f) - / (Cn(sN), f s + Xf x + nf)ds - (C(0), /) 



t 

+ ii J R f (s, e(sN))ds + o(J^j. (19) 



We are going to approximate the nonlinear term in ( fT9l) by some quantities making sense in the 
space of generalised functions since we treat the processes distributions as probability measures on 
a space D([0,T], S"(R)). Let x E C^°(R) be a non-negative function such that f R >c(y) dy = 1. 
Denote x £ (y) = x (y/e) /e, for < e < 1 and let (x e * <p(s))(x) = L x e (x — y)<p(y, s)dy be a 
convolution of a generalised function ip(s, •) with the test function >c £ (y). 

Lemma 3.2 The following uniform estimate holds 

\R f (s,e(sN)) - {{x s *( N (sN)) 2 J)\ < F 1 (e) + F 2 (eN) 
where the functions Fi and F 2 do not depend on £ and s and 

limFi(e) = lim F 2 (r) = 0. 

Proof. For definiteness we assume that x(x) = for x E (— oo, —1 — 5') U (1 + 5', oo) for some 
positive 5'. It is easy to see that if x E [k/N, (k + l)/N) for some k, then 

(j+l)/N 

W £ *( N {sN)){x) =^2^ N j(sN) J x £ {x-y)dy 

3 3/N 

= ^ y E K ^ k ~ J)/ n )Uj(sN) + 9l (N, e, x, £ N (sN)) 



where the function g%(N, e, x, £n(sN)) can be bounded as follows 

(j+i)/;v 



\ gi {N,e,x,s,Z N {sN))\ < / 



x £ (x -y)- h £ 



j/N 



< 



1 



max — 

w£[(m-l)/N,rn/N] E l 



k-j 
N 



w 

x | — 



dy 
1 

N 



o max WW 

(Ne) ^vel( m -l)/(Ne),m/(Ne)] 



2U+ £ (W\) 



Ne 
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(see the O for the notation U + ). Note that if e is fixed then U^ £ (|x'|) = 0(1) as N -> oo. 
This representation implies that 

* Cv( S iV)) 2 , /) = I £ /(Ar/JV) ^ £ - j)/iV) W*AO j + (^) 



Therefore 



e 2 (siV)) - ((x e * Cjv(siV)) 2 , /) I < J /iS (5', e, iV) + K(e, N) + O 



1 



where 



x £ ((fc - j)/N)\£ N , k (sN) - £ Nlj (sN)\ 



j:\j-k\<(l+S')eN 



and 



tf(e, AT) = Const • 



Evidently, K(e, N) = Const 



N ^ 



x. 



- 1 



and, hence, K(e, N) tends to in the limit 



"e is fixed, N -> 00". So we can include if (e, AT) into F 2 (eN). 

Consider now the term Jf tS (S', e, N). Using monotonicity of trajectories we obtain that for any 
j such that \ j — k\ < (1 + S')eN 

\^N,k( sN ) ~ £N,j(sN) \ < £ N ,k-i(l+5>)eN](sN) - £ N>k +[(l+6')eN}(sN). 

Thus we have that 

2 

N,k-[(l+5')eN]{sN) - CjV,fc+[(l+<5')eiV]( S ^))- 

k 

Integrating by parts we get the following bound 

J Ls (5', ,,A0<|£ ( fc + K 1 + ^D/AOI - |/(*. (* - [(1 + W(^) 

it 

<|E l/(*> (* + K 1 + ^D/AO - /(«, (* - [(1 + 5')eN])/N)\ 

k 

< Const ■ MD ■ (1 + S')e, 

where M = max x s \f x (s, x)\, D is a diameter of supp / x (s, x). 

Note that the last inequality is uniform in trajectories. Lemma 1X21 is proved. 



16 



A. Manita, V. Shcherbakov 



Lemma 3.3 For every 5 > 



(T) 

lim sup lim sup P^ ^ 

e^O N->oo 



R f (s, e(sN)) - ((x £ * ( N (sN)) 2 } f) ds 



>6 \ =0. 



The proof of this lemma is omitted because it is a direct consequence of Lemma 1X21 
It is easy to see that for any / 6 C^ T a map Ff t T, s (ip) '■ D([0,T], S') — > R+ defined by 



(<p(s), f s + Xf x + (if) - /i((x £ * <p(s)) 2 , /)) ds + (<p(0), f) 



is continuous, therefore for any 5 > the set {if e D([0,T], S') : Ff t T iS ((p) > 5} is open and hence 



lim sup Q)^ (<p : F f T ((p) > 5) < lim sup lim inf Q ( * { (<p : F f T e ((p) > 5) 
where is a limit point of the sequence {Qn\}n=2- Obviously that 

Qn\ (v ■ Ff,T,s(<P) >S)< P$ (sup \W' LN (t)\ > 5/2^ 



(20) 



+ AT, A 



(J2,( S , £ 2 (sA0) - ((* * Gv(siV)) 2 , /)df 



> 5/2 . 



It is easy to see that the bound © obtained in Lemma l3~T1 for the process Zf^(s) is also valid for 
the process Z'j N (s), therefore for any t 



ESi(W' f>N (t)f = <i ( J Z' LN (s)ds) < 



Ct 



since Vl N (t) is a martingale. Kolmogorov inequality implies that for any 5 > 

T 



F$ (sup |w}^(t)| > 5) < r 2 <i(w; i7V (T)) 2 = <r 2 <| I z} i7V 



CT 



(21) 



The second term in (l2"0l vanishes to zero by Lemma 1131 as iV — ► 00 and e — > 0. Therefore for any 

/ e C^ T and 5 > 

lim sup gf } : F fiT ^<p) > 5) = 0. (22) 
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Let us prove that we can replace the convolution in Ff^e by its limit which is well defined with 
respect to the measure . First of all we note that for any C > B C (R) the set of measurable 
functions h such that ||/i||oo < C is a closed subset of 5"(R) in both strong and weak topology. 
Indeed, consider a sequence of functions g n G 5c n > 1 and assume this sequence converges in 
S' to some tempered distribution G E S'. We are going to show that this generalized function is 
determinated by some measurable function bounded by the same constant C. It is easy to see that 
for every n > 1 



g n (x)f(x)dx 



<C\\f\\ L i, feS, 



where || • || L i is a norm in L 1 the space of all integrable functions. So the limit linear functional G 
on S is also continuous in L^norm 

\G(J)\<C\\f\\r>, feS. 

The space S is a linear subspace of L 1 therefore by Hahn-Banach Theorem (Theorem III-5, ifTHlO 
the linear functional G can be extended to a continuous linear functional G on L 1 with the same 
norm and such that G\s = G. Using the theorem about the general form of a continuous linear 
functional on L 1 ([ 18|) we obtain that 



G(f) = J g(x)f(x)dx, f e L\ 



where g is a measurable bounded function. Obviously that ||^||oo < C. 

Obviously that for any iV > 2 and fixed t G [0, T] we have that ((p(t, ■) G B x ) = 1, where 
(p(t, •) = if {t) is a coordinate variable on D([0, T], S'). Therefore if some subsequence {Q N ' X } °f 
{Q ( n \}n=2 converges weakly to a limit point Q X T \ then for any fixed t G [0, T] 

Q { P (ip(t) G 50 > limsupQ^ (ip(t) G B x ) = 1, (23) 

since 5i is closed. Next lemma gives an important property of the convolution on the set of 
bounded functions L°°. 

Lemma 3.4 Fix ip G L°°. Then 
1 ) for any < s < T 

((x £ *ip(s)) 2 J)- [ ip 2 (s,x)f(x)dx^0 (e-0) 



R 
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2) for any <t<T 



((x £ *^( s )) 2 ,/)- ip 2 {s,x)f{x)dx\ds^O (e-0). 



Proof of Lemma \3.4\ 



((x e * ^(s)) 2 ,/) - / p 2 (s,x)f(x)dx 



< 



< 2\\(p\ 



((x £ *<p(s))(x)\ ~(p 2 (s,x) 
((x e * y>(s))(x) - (p(s,x) 



\f(x)\ dx 
\f(x)\ dx 



To get the last inequality we used the identity a 2 — b 2 = (a + b) (a — 6) and the fact that ||x e * 
y||oo < llvlloo- To finish the proof it suffices to apply a well-known result about convergence of 
(x e * v 9 )( s ! ') to <^(s) in L^. This proves the first statement of the lemma. To get the second 
statement we use again the boundedness of and x e * ip and apply the Lebesgue theorem. Lemma 
I3.4l is proved. 

(T) 

On the set supp Q\ we can define a functional 



Ff tT (<p) = / ( V (s), f s + Xf x + fif) - ^{s), f))ds + (p(0), /). 



The equation d2"31) and Lemma IXU yield that for any <p G suppQ^ Ff,T,e(<f) ~^ F® T (ip) as 
£ — > 0. This implies that for any 5\ > 

QP {<p : \F f>T , £ (<p) - Ff >T ((p)\ >5 1 )^0 (e -> 0). 

Combining this with (l2"2l) we get that a limit point of the sequence {Q^ \}^=2 * s concentrated 
on the trajectories cp(t, •), t G [0, T], taking values in the set of regular bounded functions and 
satisfying the following integral equation 



(<p(s), f s + Xfx + ///) ~ f))ds + (p(0), /) = 



for any / E C^ T . It means that each such a trajectory is a weak solution of the equation @ in the 
sense of Definition 12. 1[ 
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3.5 Uniqueness of a weak solution 

The first order equation. Using the method of Theorem 1 in the celebrated paper ifTTl of Oleinik 
we will show here that for any measurable bounded initial function tp(x) there might be at most 
one weak solution of the equation @ in the sense of Definition 12.11 and no entropy condition is 
required. 

Let u(t, x) and v(t, x) be two weak solutions of the equation © in the region [0, T] x R with 
the same initial condition ip (not necessarily from H). Definition 12 . 1 1 implies that 

T 

r ((u(t,x) - v(t,x))(f t (t,x) + Xf x (t,x) + - u(t,x) -v(t,x))f(t,x))dxdt = 0, (24) 

R 

for any / 6 C^ T . Consider the following sequence of equations 

f t (t, x) + Xf x (t, x) + g n {t, x)f{t, x) = F{t, x), (25) 

with any infinitely differentiable function F equal to zero outside of a certain bounded region, 
lying in the half -plane t > 5i > 0, where 5% is an arbitrary small number. The functions g n (t, x) 
are uniformly bounded for all x,t,n > 1 and converges in L\ oc to the function g(t, x) = — 
u(t,x) — v(t,x)) as n — ► oo. The solution f n (t,x) G C^ T of the equation d25t is given by the 
following formula (formula (2.8) in [17]) 

f if 1 

f n (t, x) = I F(s, x + X(s — t)) exp < g n (r,x + A(r — t))dr > ds. 

T It J 

Equation (B4l yields that 

(»(«,,) - „( t , x,)F( t , I)(taft = / /( 9(i , ,) - gA .))/(«, (26) 

R + R R + R 

The right side of (126b is arbitrary small for sufficiently large n and, since the left side of d26T) does 
not depend on n, so it is equal to zero. Therefore u = v , since F is arbitrary. 

An existence of a weak solution of the equation © follows from the general theory for quasi- 
linear equations of the first order (for example, Theorem 8 in |[T7l0 . In the particular case of the 
equation © it is possible to obtain an explicit formula for a weak solution. First of all we note that 
the Cauchy problem 




Ut(t, x) = — Xu x (t, x) + fi(u 2 (t, x) — u(t, x)), u(0, x) — ijj(x), (27) 
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has a unique classical solution, if ip G C 1 (R) and there is an explicit formula for this solution. 
Indeed, using substitution u°(t, x) = u(t, x — At) we transform the equation d2"7b into the equation 

u° t {t, x) = —jiu°{t, x)(l — u°(t, x)), u°(0, x) = ip(x). 

Considering x as a parameter we obtain an ordinary differential equation which is solvable and the 
solution is given by: 

u (t, x) = J\ } — r. (28) 

So, if ip G C 1 (R), then a unique classical solution of the equation d2"7b is given by the following 
formula 

X) ~ l-V(x + At)+^(a; + At) e -M- (29) 

If we approximate any measurable bounded function g in L\ oc by a sequence of smooth func- 
tions {g n , n > 1}, then the sequence of corresponding weak solutions {u n (t, x),n > 1}, where 
u n (t, x) is defined by the formula with ip = g n , converges in L\ oc to the weak solution of the 
equation with initial condition g by Theorem 11 in ifTTl or Theorem 1 in ifTOll . It is easy to show 
by direct calculation that the L^-limit of the sequence {u n (t, x),n > 1} is given by the same 
formula (l2"9t with tp = g. 

The formula d2"9t yields that if ip e H(R), then u(t, •) G H(R) as a function of x for any 
fixed t > 0. If a function u(t,x) is a weak solution of the equation (|2"71) . then this function is 
differentiable at a point (i, x) iff the initial condition ip(y) is differentiable at the point y = x + Xt. 

KPP-equation. The equation © is a quasilinear parabolic equations of the second order and is 
a particular case of the famous KPP-equation in [9|. It is known that there exists a unique weak 
solution u(t, x) of the problem 

Ut(t, x) = ^u xx (t, x) + n(u 2 (t, x) — u(t, x)), u(0, x) = ip(x), 

for any bounded measurable initial function ip, this solution is in fact a unique classical solution 
and if ip G H (R), then u(t, •) G H(R) as a function of x for any fixed t. We refer to the paper ifTTl 
for more details. 



4 System with fixed number of particles 

In this section we deal with the situation "N is fixed, t — > oo". 
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4.1 Proof of Theorem El 

Let a(y,w) be the rate of transition from the state y = (yi, . . . , y^) G T to the state w = 
(wi, . . . ,w N ) G T for the Markov y(t) chain. Define a(y) = ° r (y^ w )- From definition of 

the particle system it follows that 



N 



{yi>yj}- 



Since ^{yi>yj} — N(N — l)/2 we have uniformly in y G T 

(m) 

cr*,N < <?(y) < <7*n (30) 

with a* t N = (a + (3)N and cr^ = (a + /3)JV + ^n{N - l)/2. A discrete time Markov chain 
n = 0, 1, . . .} on the state space S with transition probabilities 

— y =fc yj 

°{y) ' (3D 
0, y = w, 

is an embedded Markov chain of the continuous time Markov chain (y(t), t > 0). 
Theorem 12 .21 is a consequence of the following statement. 

Lemma 4.1 The Markov chain {Y(n), n = 0, 1, . . .} is irreducible, aperiodic and satisfies to the 
Doeblin condition: there exist e > 0, itlq G N and finite set A G T such that 

P{Y(m )eA\Y(0) = Y } > e, (32) 

for any Y G S. Therefore this Markov chain is ergodic ([5]). 

Proof of Lemma W7\ We are going to show that condition d3*2l holds with A = {(0, . . . ,0)}, 
m = N, 



min (a, (3, [xn/N) \ N 
°*n J 



> 0. 



'N 

The transition probabilities of the Markov chain {Y(ri), n = 0, 1, . . .} are uniformly bounded from 
below in the following sense: if a pair of states (z, v) is such that a(z, v ) > (or, equivalently, 
p(z, v) > 0) then d3lT) implies that 

p(z, v) > min (a,/3,/i N /N) /a* N . (33) 
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So to prove (l3~2l we need only to show that for any y there exists a sequence of states 

v° = y, v\ v 2 , v N = (0,...,0) (34) 

which can be subsequently visited by the Markov chain {Y(n), n = 0, 1, . . .}. The last means that 
i.e. p(u n-1 , v n ) > for every n = 1, . . . , N and hence 

/min(a,/3, /iAr/iV) x a 



P {Y(N) G A | Y(0) = y} > Y[p(v n -\v n ) > 



n=l V °*N 

as a consequence of the uniform bound (l33l) . 

To prove existence of the sequence 434b let us assume first that y = (yi, . . . , y^) ^ 0. Choose 
and fix some r such that y r = maxj/i > 0. Denote by 

i 

n = # {j : yj = 0} 

the number of left-most particles. Let the right-most particle y r move n steps to right: 

v n _ v n-i = e (N) ^ n =l,...,no. 

This can be done by using of jumps to the nearest right state. So Y(n ) = v n ° has exactly n 
particles at and N — n particles out of 0. Denote by i ri0 +i < K +2 < ••• < indices of 
particles with t>™ > 0, a = no + 1, . . . , N. Let now the Markov chain Y transfer each of these 
particles to 0: 

v a -v a - x = -v^e { "\ a = n + l,...,N. 



It is possible due to transitions provided by the interaction. 

To complete the proof we need to consider the case y = [y Xl . . . , y^) = 0. It is quite easy: 

^ = eS N) , v 2 = 2e[ N \ v 3 = 3e[ N \ v N ~' = (N - l)e[ N \ v N = 0. 

Proof of the lemma is over. 

Denote by 7r y = (7r y (y), y = . . . ,y N ) G T) a unique stationary distribution of the Markov 
chain {Y(n), n = 0, 1, . . .}. 

The proof of Theorem \2.2\ is now easy. First of all let us show that the uniform bound a(y) > 
a^N implies existence of a stationary distribution for the Markov chain y(t). Indeed, it is easy 
to check that if n Y is the stationary distribution of the embedded Markov chain Y and Q is the 
infinitesimal matrix for the chain y(t), then a vector with positive components s = (s(w),w G S) 
defined as 

7T Y (w) 

s H = / \ , 

a(w) 
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satisfies to the equation sQ = 0. So for existence of a stationary distribution of the chain y(t) it is 
sufficient to show s(w) < +00. It is easy to check the last condition: 



w€S 



^ ^ 71 (W) 1 ^ Y 1 
> S{W) = > r -^~ < > TT (W) = — . 

Therefore the continuous-time Markov chain (y(t), t > 0) has a stationary distribution n = (n(y), y G T) 
of the following form 

Ay) = 



E 



71 (W) 



a(w) 

Denote p yw (t) = P {y(t) = w \ y(0) = 0}. The next step is to prove that the continuous-time 
Markov chain y{t) is ergodic. To do this we show that the following Doeblin condition holds: for 
some jo G S there exists h > andO < 5 < 1 such thatpij (h) > 5 for alii G S. It is well-known 
((5 ]) that this condition implies ergodicity and moreover 



\ Pij (t)-7r(y)\ < (1-5) 



[t/h] 



Let Tk, k > 0, be the time of stay of the Markov chain y(t) in fc— th consecutive state. Condition 
on the sequence of the chain states y k , k > 0, the joint distribution of the random variables r k , k > 
0, coincides with the joint distribution of independent random variables exponentially distributed 
with parameters a(y k ), k — 0, 1, . . . , n, so the transition probabilities of the chain y(t) are 

„ n 

PyUt) = EE P {^ ffi )W e -H(*-*») ft aiy^e-^-^-^ dt, . . . dt n , 

n (y-^w) k=l 

where n corresponds to the number of jumps of the chain y during the time interval [0, t], the inner 
sum is taken over all trajectories (y — > w) = {y = y , yi, . . . , y n = w} with n jumps, integration 
is taken over = {0 = t < t\ < ■ ■ ■ < t n < t}, and 

P {(y -> w)} = p(y, yi)p(yi, y 2 ) ■ ■ -p(yn-i, w) 

is a probability of the corresponding path for the embedded chain. The equation (l3H implies that 
the integrand in p yw (t) is uniformly bounded from below by the expression 

((j, i 7v) n exp(-cr^t 1 ) • • • exp(-a* N (t n - t n - X )) exp(-a* N (t - t n )), 
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and, hence, 



/n 
e -a{w){t-t n ) TT a^Ae-^Vk-Mk-tk-i) ^ . . . dt n > 

An n- — -L 



P{U t = n} e-i^-^y, 



where H t is a Poisson process with parameter a* t N- It provides us with a lower bound for the 
transition probabilities of the time-continuous chain: 



Pyz 



It is easy now to get a lower bound for probabilities p n (y,z ). Fix some z E S and denote 
£ = 7i Y (z ). It follows from ergodicity of the chain F, that for any fixed zq 

p m (y,z )>n Y (z )/2 = Z/2>0, 

for all m>mi= mi(y). For a Doeblin Markov chain we have more strong conclusion, namely, 
the above number mi does not depend on y. Let us fix such m x and show that the continuous-time 
Markov chain y(t) satisfies to the Doeblin condition. Indeed, 

Pyzo (t) > ( J2p n (y^o)P{^t = n}+ J2p n (y^o)P{^t^n})e-(<-^y 

\n<mi n>m\ / 

> p"(y^o)P{n t = n} e -0*- ff '."> 

n>mi 

2 

Hence, the Doeblin condition holds: we choose any z as j , take a corresponding m x , fix any 
/i > and put 

§ = i p {n h > m 1 }e-(' T "- a *' N ) h . 

Proof of the theorem is over. 



4.2 Evolution of the center of mass 

Consider the following function on the state space Z N : m(xi, . . . , xn) = (xi + • • • + xn) /N. So 
if each particle has the mass 1 and x±(t) : . . . , xjv(t) are positions of particles, then m(x 1 (t), 
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is the center of mass of the system. We are interested in evolution of E m(x ly t), . . . , x N (t)). A 
direct calculation gives that 



[<x - P)-^rk^2\ x i- x A, (35) 



N 2 

Kj 



where we have used the following equalities 



m (x ± e[ N ^ — m (x) = i"r^> 
;;/ i x — (xi — Xj)e\ N ^ — m{x) = 1 — J 



Note that the summand 



N 2 

Kj 



X j\ 



added by the interaction to the "free dynamics" drift (a—/3) depends only on the relative disposition 
of particles. So the center of mass of the system moves with speed which tends to the value 

N -1 

(a - (3) - fj, N E vr Fi - x 2 \ 

as t goes to infinity. Here \x\ — x 2 | is the mean distance between two particles calculated with 
respect to the stationary measure ir of the Markov chain Y. 

Using this fact and Theorem l2.2l we can describe the long time behavior of the particle system in 
the initial coordinates x's as follows. Theorem 12 . 21 means that the system of stochastic interacting 
particles possesses some relative stability. In coordinates y the system approaches exponentially 
fast its equilibrium state. In the meantime the particle system considered as a "single body" moves 
with an asymptotically constant speed. The speed differs from the mean drift of the free particle 
motion and this difference is due to the interaction between the particles. 



5 On travelling waves and long time evolution of solution of 
PDE 



We deal here with partial differential equation in variables (t,x) G R + x R. 
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Definition 5.1 Function w = w(x) is called a travelling wave solution of some PDE if there exists 
v G R such that the function u(t, x) = w(x — vt) is a solution of this PDE. The number v is speed 
of the wave w. 

We are interested only in the travelling waves having the following properties: Ul) w(x) G [0, 1]; 
U2) w(x) and dw(x)/dx have limits as x — > ±oo and, besides, w(— oo) = 1 and w(+oo) = 0. We 
identify two travelling waves w\(x) and W2(x) if Wi(x) = ^(x — c) for some c. 

For any probabilistic solution < u(t, x) < 1 we define a function r(t) such that u(t, r(t)) = 
|. Let a function be a travelling wave solution. Without loss of generality we can assume 
thatw(O) = 1/2. 

Definition 5.2 A solution u(t, x) converges in form to the travelling wave w(x) if 

lim u(t, x + r(t)) = w(x), 

uniformly on any finite interval. The solution u(t,x) converges in speed to the travelling wave 
w(x) if there exists r'(t) = dr(t) / dt and lim r'(t) = v, where v is a speed of the travelling wave 

t— >+oo 

w(x). 

First order equation. It is easy to check for the equation © for every v < X there exists a unique 
(up to shift) travelling wave solution having properties U1-U2 and this travelling wave solution is 

given by the following formula wJx) = 1 1 + exp ( — — — x 

Proposition 5.1 If for some C > 0, v > 0, the initial profile ift(x) G H(R) of the equation (0) has 
the following asymptotic behavior 1 — ip(x) ~ C exp(ux), as x — > — oo, then there exists xq G R 
such that for every x G R 

\u(t, X) — W v (x — Xq — vt) I — > 0, 

as t — >• oo, where v — X — p,jv. 

The proof of Proposition l5. H is a direct calculation based on the exact formula (l29l) . 

The formula (l28t yields that u°(ti,x) > u°(t 2 , x) for any t\ < t 2 and this observation imme- 
diately implies the following statement: for every x u°(t, x) — > I{y.My)=i}(x), as t — > oo. As a 
direct application of this property we can obtain the following 

Proposition 5.2 Assume that the initial profile has the form i/j(x) = I{ y< b}(x) for some b G R. 
Then 

\u(t, x) — I{ y< b}(x — At) I — > 0, t — > oo. 

So the function w(x) = I{ y < y(x) is a unique (up to shift) non-increasing continuous from the 
right travelling wave corresponding to the maximal possible speed v = X. It is easy to see that this 
function is a limiting case of w v (x) as v — > A — 0. 
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Second order equation. The existence of travelling waves for parabolic partial differential equa- 
tions was a subject of studying in many papers followed to the paper [9|. A review of many results 
can be found in ll2"4l (see also (20 1) and for completeness of the text we mention some of them. 
Reformulating the well-known results (|24|) we obtain that travelling waves of the equation © 
can move only from the right to the left. It means that the speed of any travelling wave is negative 
and, moreover, is bounded away from 0. 

Proposition 5.3 For equation every v < v * = — ^/47/i there exist and unique (up to shift) 

travelling wave solution with speed v. There are no other travelling wave solutions satisfying the 
conditions Ul and U2. 

If a function / = f(x) is such that f(x) < 1, f(x) — > 1 as x — ► —00 and there exists a limit 
x = lim x _1 log(l — f(x)) > 0, then the number x is called Lyapunov exponent of the function 

X— >— 00 

f (at minus infinity). It is well known (EU) that for the equation © a travelling w(x) with speed v 
has the following Lyapunov exponent at minus infinity 

x(t>) = ( — v — \/ v 2 — 47/i) / (27) . 



Hence we get that for the travelling wave with minimal in absolute value speed = — \f%rfji the 
Lyapunov exponent is x(i>*) = a/ Jtpf- 

Proposition 5.4 (|24|) Assume that an initial function ip(x) has a Lyapunov exponent x. Then 

a) if *c > a//V7 t nen tne solution u(t, x) of the problem converges inform and in speed to 
the travelling wave moving with the minimal speed v* = —y/Ajp; 

b) if x < a//U/7 then the solution u(t, x) of the problem (0) converges inform and in speed to 

the travelling wave with speed v = — yjH + —j , or, in other words, x(t>) = x. 

We see from the above analysis that both the first order PDE © and the second order PDE © 
exhibit similar long-time behavior of their solutions. This seems very natural if we recall from The- 
orem EH] that the both equations arise as hydrodynamical approximations of the same stochastic 
particle system. 

A Appendix 

A.l Strong topology on the Skorokhod space 

Remind that Schwartz space S'(R) is a Frechet space (see ifTHI ). In the dual space S"(R) of tem- 
pered distributions there are at least two ways to define topology (both not metrizable): 
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1) weak topology on 5"(R), where all functionals ( • , 0) , <fi G S(R) are continuous. 

2) strong topology (s.t.) on 5"(R), which is generated by the set of seminorms 

p A (M) = sup \{M,<f>)\ : A C S(R) - bounded 

Below we shall consider S'(R) as equipped with the strong topology. The problem of introducing 
of the Skorokhod topology on the space D T (S') := D([0, T], S'(R)) was studied in [15] and Q. 
We follow these papers. For each seminorm p A on S'(R) we define the following pseudometric on 

D([0,T],S'(R)) 



d A (y, z) = inf <^ sup p A (y t - z x{t) ) + sup 



log 



X(t) - X(s) 



y,zeD([0,T],S'(R)), 



where inf is taken over the set A = {A = A(t), t G [0, T}} of all strictly increasing continuous 
mappings of [0, T] onto itself. Introducing on D([0, T], 5"(R)) the projective limit topology of 
{d A (-, •)} we S et a completely regular topological space. 



A.2 Mitoma theorem 

Let Bd t (s') be the corresponding Borel a-algebra. Let {P n } be a sequence of probability measures 
on (D t (S'),B Dt ( S >))- For each G 5(R) consider a map : y G D T (S') (y,(j>) G D T (R). 
The following result belongs to I. Mitoma [fT51. 

Theorem A.l Suppose that for any G 5(1?) sequence {PnZ^ 1 } w m Dt(R). Then the 
sequence {P n } itself is tight in Dt(S'). 

A.3 Probability measures on the Skorokhod space: tightness 

Let G [0,T])} ngN be a sequence of real random processes which trajectories are right- 

continuous and admit left-hand limits for every < t < T . We will consider £ n as random 
elements with values in the Skorokhod space D T (R) := D ([O^R 1 ) with the standard topology. 
Denote Pj. the distribution of defined on the measurable space (D T (R), B (D T (R))). The 
following result can be found in ifD . 

Theorem A.2 The sequence of probability measures {P£ } ngN is tight iff the following two con- 
ditions hold: 

l)for any e > there is C(e) > such that 

supP? ( sup \C\>C(e)) <e; 

n \0<t<T / 
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2) for any e > 

lim lim sup P% (£. : w'{£\ 7) > e) 



where for any function f : [0, T] — > R and any 7 > we define 



7) = inf max sup \f(t) 

{*i};=i J<r i l < s <t<t i+1 



moreover the inf is over a// partitions of the interval [0, T] such 



that 



= t Q <ti<---<t r = T, 



U ~ U-i > 1, i = l,...,r. 



The following theorem is known as the sufficient condition of Aldous [8, Proposition 1.6]. 
Theorem A.3 Condition 2) of the previous theorem follows from the following condition 



where TZt is the set of Markov moments (stopping times) not exceeding T . 
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